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We present evidence for the saturation of the Froissart bound at high energy for all hadronic total 
cross sections at high energies, and use this to unify pp (and pp) total cross sections over the energy 
range from cyclotrons to colliders to ultra-high energy cosmic rays, an energy span from \/s = 4 
GeVto80TeV. 



Introduction. High energy cross sections for the scat- 
tering of hadrons should be bounded by a ~ In 2 s, where 
s is the square of the cms energy. This fundamental re- 
sult is derived from unitarity and analyticity by Froissart 
I], who states: 

"At forward or backward angles, the modulus of the 
amplitude behaves at most like s In 2 s, as s goes to infin- 
ity. We can use the optical theorem to derive that the 
total cross sections behave at most like In 2 s, as s goes to 
infinity" . 

In this context, saturating the Froissart bound refers 
to an energy dependence of the total cross section rising 
no more rapidly than In 2 s. 

It will be shown that the Froissart bound is saturated 
at high energies in ■yp, ir^p and Jpp and pp scattering 
[H, as well as in j*p scattering Q, as seen from Deep 
Inelastic Scattering (DIS) in e + p — > e + X. 

Using Finite Energy Sum Rules (FESR) derived from 
analyticity constraints — in order to anchor accurately 
cross sections at cyclotron energies [3] — we will make 
precise predictions about the total pp cross section, the 
p-value (the ratio of the real to the imaginary portion 
of the forward pp scattering amplitude), as well as the 
shape of the differential elastic scattering cross section, 
durmeei/dt, at the LHC. Further, we will make predic- 
tions of the total pp cross section at cosmic ray energies, 
up to 50 TeV, and will compare them to the latest ex- 
periments. 

Data selection. We make the following major assump- 
tions about the experimental data that we fit: 

1. The experimental data can be fitted by a model 
which successfully describes the data. 

2. The signal data are Gaussianly distributed, with 
Gaussian errors. 

3. The noise data consists only of points "far away" 
from the true signal, i.e., "outliers" only. 

4. The outliers do not completely swamp the signal 
data. 



We will use the "Sieve" algorithm [5| to remove "out- 
liers" in the cross section and p- values that we will fit, in 
order to improve the accuracy of our fits. The "Sieve" 
Algorithm does the following: 

1. Make a robust fit of all of the data (presumed 
outliers and all) by minimizing Aq, the Lorentzian 
squared, defined as 

N 

A 2 (a; x) = ]T In {l + 0.179A X 2 (x i; a)} , 

i=l 

where A X 2 (^;a) ^ ( ^1^1^ . (1) 

The M -dimensional parameter space of the fit is 
given by a = {ai,.. . ,a A /}; x = {x 1} .. .,x N } is 
the abscissa of the N experimental measurements 
y = {yi,...,y N } that are being fit; y(x 2 ;a) is 
the theoretical value at Xi and 0$ is the experi- 
mental error. Minimizing Aq gives the same total 
Xmin = Eili Ax?(x<; a) as that found in a x 2 fit, 
as well as rms widths (errors) for the parameters — 
for Gaussianly distributed data — that are almost 
the same as those found in a x 2 fit- The quanti- 
tative measure of "far away" from the true signal, 
i.e., point i is an outlier, is the magnitude of its 

A x f(^;«) = ( ^f^ ) 2 . 

If Xmi n is satisfactory, make a conventional x 2 fit 
to get the errors and you are finished. If Xmin * s 
not satisfactory, proceed to step [2] 

2. Using the above robust Aq fit as the initial estima- 
tor for the theoretical curve, evaluate Ax 2 (xj;a), 
for each of the N experimental points. 

3. A largest cut, A X f (xj; a) max , must now be se- 
lected. For example, we might start the process 
with A X f (xf, «) max = 9. If any of the points 
have AXi( x H a ) > ( X H «)max, reject them— 
they fell through the "Sieve". The choice of 
^Xi( x i> a) max is an attempt to pick the largest 



2 



"Sieve" size (largest Axf (£»; a)max) that rejects all 
of the outliers, while minimizing the number of sig- 
nal points rejected. 

4. Next, make a conventional x 2 fit to the sifted set — 
these data points are the ones that have been re- 
tained in the "Sieve". This fit is used to esti- 



mate xi 



Since the data set has been trun- 



cated by eliminating the points with Axf (aij; a) > 
Axf {xi\ a)max, we must slightly renormalize the 
X min found to account for this, by the factor 
■£=1.027, 1.14, 1.291 for Axf max = 9,6,4. If the 
renormalized Xmim *- e -> ^ x Xmin i s acceptable — 
in the conventional sense, using the x 2 distribution 
probability function — we consider the fit of the data 
to the model to be satisfactory and proceed to the 
next step. If the renormalized Xmin i s n °t accept- 
able and Ax 2 (x{; a) max is not too small, we pick a 
smaller Axf (a:*; a)max and go back to step El The 
smallest value of Axf (a:*; ct) max that makes much 
sense, in our opinion, is Axf (a;*; a) max > 2. One of 
our primary assumptions is that the noise doesn't 
swamp the signal. If it does, then we must dis- 
card the model — we can do nothing further with 
this model and data set! 

5. From the x 2 fit that was made to the "sifted" data 
in the preceding step, evaluate the parameters a. 
Next, evaluate the M x M covariance (squared er- 
ror) matrix of the parameter space found in the x 2 
fit. We find the renormalized squared error ma- 
trix of our x 2 fit by multiplying the covariance ma- 
trix by the square of the factor r x 2 (we find r x 2 ~ 
1.02,1.05, 1.11 and 1.14 for Axf (a*; a) max = 9, 6, 
4 and 2). The values of r x 2 > 1 reflect the fact 
that a x 2 fit to the truncated Gaussian distribution 
that we obtain has a rms (root mean square) width 
which is somewhat greater than the rms width of 
the x 2 fit to the same untruncated distribution. Ex- 
tensive computer simulations demonstrate that this 
robust method of error estimation yields accurate 
error estimates and error correlations, even in the 
presence of large backgrounds. 

You are now finished. The initial robust A 2 , fit has been 
used to allow the phenomenologist to find a sifted data 
set. The subsequent application of a x 2 fit to the sifted 
set gives stable estimates of the model parameters a, as 
well as a goodness-of-fit of the data to the model when 
Xmin i s renormalized for the effect of truncation due to 
the cut Axf {xi\ a) max . Model parameter errors are found 
when the covariance (squared error) matrix of the x 2 fit 
is multiplied by the appropriate factor (fV 8 ) for the cut 
Ax 2 (x l ;a) max . 

Analyticity constraints on hadronic cross sections. 
Block and Cahn Q used an even real analytic amplitude 
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where v is the nucleon (pion) laboratory energy and m 
is the proton mass. Using the optical theorem, the even 
cross section is 

a+(v) = c + ci ln(p/m) + c 2 ln 2 (^/m) + j3-pi{v /m) M_1 , 

(4) 

where here the coefficients Co, c±, c 2 and /3p/ have di- 
mensions of mb. Since — = for v ^> m, i.e., for 
the high energy regime in which we are working, we will 
henceforth refer to equations similar to that of Eq. (Q} as 
In 2 s fits if c 2 ^ and as Ins fits if c 2 = 0. 

We now introduce f+(v), the true even forward scatter- 
ing amplitude (which of course, we do not know!), valid 
for all i/, where f+{v) = [fpp(v)+ fpp(v)]/2, using forward 
scattering amplitudes for pp and pp collisions. Using the 
optical theorem, the imaginary portion of /+ {v) is related 
to the even total cross section o~ even (v) by 



47T 



for v > m. 



(5) 



Next, define the odd amplitude vf + {v) as the differ- 
ence 



v f+(y) = v f+iy) - /+M > 



(6) 



which satisfies the unsubtracted odd amplitude disper- 
sion relation 

TT I _ it I & 1/^. 



Since for large v, the odd amplitude vf + {y) ~ v a {a < 0) 
by design, it also satisfies the super-convergence relation 



Im vf+(v) dv = 0. 



In Ref. [8|, the FESRs are given by 



(8) 



v n Im/ dv = . 



a+n+1 



a + n + 1 



, n = 0, 1, . . . , oo, 



(9) 

where f(v) is crossing-even for odd integer n and 
crossing-odd for even integer n. In analogy to the n = 1 
FESR of Ref.. [8], which requires the odd amplitude 
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vf{v), Igi and Ishida inserted the super-convergent am- 
plitude of Eq. (|6]) into the super-convergent dispersion 
relation of Eq. flSJ), obtaining 



v Im 



f + (u)-f + (u) 



dv. 



(10) 



We note that the odd difference amplitude vim f+(v) sat- 
isfies Eq. ([HJ , a super-convergent dispersion relation, even 
if neither vim f+{v) nor vim f+(v) satisfies it. Since 

the integrand of Eq. (fTU|) . ^Im f + (v) — f+(v) , is super- 
convergent, we can truncate the upper limit of the inte- 
gration at the finite energy vq, an energy high enough 
for resonance behavior to vanish and where the differ- 
ence between the two amplitudes — the true amplitude 
f+(v) minus /+(;/), the amplitude which parametrizes 
the high energy behavior — becomes negligible, so that 
the integrand can be neglected for energies greater than 
vq. Thus, after some rearrangement, we get the even 
finite energy sum rule (FESR) 



vim f+{v) dv 



vim f+(v)dv. (11) 



Next, the left-hand integral of Eq. (fTTj) is broken up 
into two parts, an integral from to m (the 'unphysi- 
cal' region) and the integral from m to vq, the physical 
region. We use the optical theorem to evaluate the left- 
hand integrand for v > m. After noting that the imagi- 
nary portion of f+{v) — for < v < m, we again use 
the optical theorem to evaluate the right-hand integrand, 
finally obtaining the finite energy sum rule FESR(2) of 
Igi and Ishida [6], in the form: 



I vlmf + (v) dv + — vp a cvcn (v) dv = 
Jo 4?r J m 



-in 



vpa+(v) dv. 



(12) 



We now enlarge on the consequences of Eq. ([12")) . We 
note that if Eq. (|12[) is valid at the upper limit vq, it 
certainly is also valid at vq + Avq, where Avq is very 
small compared to vq, i.e., < Avq -C vq. Evaluating 
Eq. (|12l) at the energy vq + Avq and then subtracting 
Eq. (fT2)) evaluated at vq, we find 
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vpa+(v)dv. (13) 



Clearly, in the limit of Avo — !• 0, Eq. (113)) goes into 

o-cvon(^o) = v + (v ). (14) 
Obviously, Eq. (TH1) also implies that 

0"cvcn(^) = G+{v) for all V > Vq, (15) 



but is most useful in practice when vq is as low as possi- 
ble. The utility of Eq. (1151) becomes evident when we rec- 
ognize that the left-hand side of it can be evaluated using 
the very accurate low energy experimental crossing-even 
total cross section data, whereas the right-hand side can 
use the phenomenologist's parameterization of the high 
energy cross section. For example, we could use the cross 
section parameterization of Eq. (U) on the right-hand side 
of Eq. (IT51) and write the constraint 



[crpp(v) + crppiy)} /2 = c + c\ miyjra) + c 2 m 2 (v/m) 



+P v ,(v/m) 



(16) 



where a pp and ap p (v) are the experimental pp and pp 
cross sections at the laboratory energy v. Equation (TTi)) 
(or Eq. (|15p) is our first important extension, giving us an 
analyticity constraint, a consistency condition that the 
even high energy (asymptotic) amplitude must satisfy. 

Reiterating, Eq. (|15|) is a consistency condition im- 
posed by analyticity that states that we must fix the even 
high energy cross section evaluated at energy v > v (us- 
ing the asymptotic even amplitude) to the low energy 
experimental even cross section at the same energy v, 
where Vq is an energy just above the resonances. Clearly, 
Eq. (frjj also implies that all derivatives of the total cross 
sections match, as well as the cross sections themselves, 
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(v), n = 0,l,2,... v> i/ 0) (17) 



giving new even amplitude analyticity constraints. Of 
course, the evaluation of Eq. (TTT)) for n — and n = 
1 is effectively the same as evaluating Eq. (TT7)) for n = 
at two nearby values, vq and v\ > vq. It is up to 
the phenomenologist to decide which experimental set of 
quantities it is easier to evaluate. 

We emphasize that these consistency constraints are 
the consequences of imposing analyticity, implying sev- 
eral important conditions: 

1. The new constraints that are derived here tie to- 
gether both the even hh and hh experimental cross 
sections and their derivatives to the even high en- 
ergy approximation that is used to fit data at en- 
ergies well above the resonance region. Analyticity 
then requires that there should be a good fit to the 
high energy data after using these constraints, i.e., 
the x 2 per degree of freedom of the constrained fit 
should be ~ 1, if the high energy asymptotic am- 
plitude is a good approximation to the high energy 
data. This is our consistency condition demanded 
by analyticity. If, on the other hand, the high 
energy asymptotic amplitude would have given a 
somewhat poorer fit to the data when not using the 
new constraints, the effect is tremendously magni- 
fied by utilizing these new constraints, yielding a 
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very large \ 2 P er degree of freedom. As an exam- 
ple, both Block and Halzen Q and Igi and Ishida 
p conclusively rule out a In s fit to both n^p and 
pp and pp cross sections and p- values because it has 
a huge x 2 per degree of freedom. 

2. Consistency with analyticity requires that the re- 
sults be valid for all v > i>q, so that the constraint 
doesn't depend on the particular choice of v. 

3. The non-physical integral JJ™ flm f+{v) dv is not 
needed for our new constraints. Thus, the value of 
non-physical integrals, even if very large, does not 
affect our new constraints. 

Having restricted ourselves so far to even amplitudes, 
let us now consider odd amplitudes. It is straightforward 
to show for odd amplitudes that FESR(odd) implies that 



'M = ~dU*~M' u = ' 1 ' 2 '- 



dv n 



v>v , (18) 



where <r_ (y) is the odd (under crossing) high energy cross 
section approximation and cr dd(j / ) is the experimental 
odd cross section. 

Thus, we have now derived new analyticity constraints 
for both even and odd cross sections, allowing us to con- 
strain both hh and hh scattering. Block and Halzen [2j 
expanded upon these ideas, using linear combinations of 
cross sections and derivatives to anchor both even and 
odd cross sections. A total of 4 constraints, 2 even and 
2 odd constraints, were used by them in their successful 
In 2 s fit to pp and pp cross sections and p- values, where 
they first did a local fit to pp and pp cross sections and 
their slopes in the neighborhood of vq — 7.59 GeV (cor- 
responding to y/s Q = 4 GeV), to determine the experi- 
mental cross sections and their first derivatives at which 
they anchored their fit. The data they used in the high 
energy fit were pp and pp cross sections and p- values with 
energies y 7 ? > 6 GeV. Introducing the even cross section 
o~q(v), they parameterized the high energy cross sections 
and p values [IJ as 



00 



(v) = c + ci In (-) + c 2 In 2 (-) + fo, (-)" \ (19) 
\m/ Vm/ \77i/ 

a ± (u)=a (-)± S(-Y~\ (20) 
v m / \mJ 



1 [ 7T 



"Ci + C 2 



47T „ , , „ / ira \ ( v 
— f + ± J tan — - 
v \ 2 / Vra 



(21) 



We note that the even coefficients co,ci,C2 and j3-p> are 
the same as those used in Eq. (fl6l) . The real constant 
/+(0) is the subtraction constant^, [t| required at v = 
for a singly-subtracted dispersion relation. They also 
used p = 0.5. The odd cross section in Eq. (|2"0")) is given 



^(m) Q ' described by two parameters, the coeffi- 
cient 5 and the Regge power a < 1, so that the difference 
cross section between pp and pp vanishes at high energies. 

We now have new analyticity constraints for both even 
and odd amplitudes. The fits are anchored by the exper- 
imental cross section data near the transition energy vq. 
These consistency constraints are due to the application 
of analyticity to finite energy integrals — the analog of 
analyticity giving rise to traditional dispersion relations 
when it is applied to integrals with infinite upper limits. 

7P, ^p, pp and pp scattering. Using the 4 constraints 
of the previous Section, we show in Fig. [T] both the In 2 s 
fit and the In s fit of Eq. {jT9j) and Eq. ([20]) to the exper- 
imental 7r ± p total cross sections. The high energy fit is 
anchored to the very accurate data a at yfs = 2.6 GeV. 
Because of the 4 constraints, the In 2 s fit needs only to 
fit the two parameters c\ and C2, whereas the Ins fit sets 
C2 = and fits only c\. We see that the In 2 s fit, which 
saturates the Froissart bound, gives an excellent fit to the 
data, while the Ins fit is ruled out. Shown in Fig. [2] are 
the p-values for n ± p scattering: again, we see that the 
In s fit is ruled out, whereas we get a good fit when we 
saturate the Froissart bound. 

In Fig. [3] we compare the results of our fitted irp 
cross section o"o (from Eq. (|19[0 with a rescaled version 
of cr (7p) that was obtained from a fit to all known high 
energy ■yp cross sections. From about 2 < ^s < 300 
GeV, the two saturated In 2 s fits are virtually indistin- 
guishable. Thus we conclude that high energy yp total 
cross sections also go as In 2 s. 




^ s, in GeV 



FIG. 1: The fitted total cross sections cr^+ p and cF n - p in mb, vs. 
yfs, in GeV, using 4 constraints. The circles are the sieved data 
for tt~p scattering and the squares are the sieved data for 7r+p 
scattering for y/s > 6 GeV. The dash-dotted curve (tt+p) and the 
solid curve (n~ p) are In 2 s fits. The short dashed curve (tt+p) and 
the long dashed curve (ir~ p) are Ins fits and clearly do not fit the 
data. 

We next consider high energy pp and pp scattering. 
Shown in Fig. @]are fits to the total cross sections, where 
again, for the In 2 s fit, only 2 parameters, c± and c 2 are 
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FIG. 2: The fitted p-values, p„.+ p and p^—p, vs. y/s, in GeV, using 
4 constraints. The circles are the sieved data for ir~p scattering 
and the squares are the sieved data for ir + p scattering for yTs > 6 
GeV. The dash-dotted curve (ir + p) and the solid curve (tt~ p) are 
In 2 s fits. The short dashed curve (ir + p) and the long dashed curve 
(n~~p) are Ins fits, which again are bad fits to the experimental 
data. 



FIG. 4: The fitted total cross sections u pp and a pp in mb, vs. \/s, 
in GeV, using 4 constraints. The circles are the sieved data for pp 
scattering and the squares are the sieved data for pp scattering for 
y^s > 6 GeV. The dash-dotted curve (pp) and the solid curve (pp) 
are In" fits to the high energy data. The short dashed curve (pp) 
and the long dashed curve (pp) are In s fits to the high energy data. 
Clearly, the In s fits do not fit the experimental data. 
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FIG. 3: The circles are the cross section data for tt~ p scattering 
and the squares are the cross section data for 7r+p scattering, in 
mb, vs. yfs, in GeV, for all of the known data. The dashed curve 
is the In 2 s irp fit to the high energy cross section data of the even 
amplitude cross section, <ro. The solid curve is the fit of the ■yp cross 
section data , for cms energies y/s > 2.01 GeV, whereas the irp data 
(cross sections and p-values) were fit for cms energies y's > 6 GeV. 
The two fitted curves are virtually indistinguishable in the energy 
region 2 < ^s < 300 GeV. 



fit, after anchoring the high energy data to the accurate 
cyclotron data at y/s — 4 GeV. Once again, we see that 
we have an excellent In 2 s fit to all of the high energy 
total cross section data when we anchor our fit (of only 
2 parameters!) to the low energy data, whereas the Ins 
fit to these data fails. In Fig. [5]we compare the nucleon- 
nucleon p-values to our 3 parameter fit of Co, C\ and the 
subtraction constant /+(0). 

Before employing the "Sieve" algorithm on the totality 
of 212 high energy pp and pp cross sections, the x 2 /d-f- 
was 5.7 for the In 2 s fit, clearly an unacceptably high 



FIG. 5: The fitted p-values, p pp and pp p , vs. y/s, in GeV, using 
4 constraints. The circles are the sieved data for pp scattering and 
the squares are the sieved data for pp scattering for yfs > 6 GeV. 
The dash-dotted curve (pp) and the solid curve (pp) are In 2 s fits 
The short dashed curve (pp) and the long dashed curve (pp) are 
Ins fits. Clearly, the Ins fits do not fit the experimental data. 



value. After "sieving", the renormalized x 2 /rf./. was 
1.09, for 184 degrees of freedom, using a A% 2 > 6 cut. 
The total \ 2 was 201.4, corresponding to a probability of 
fit 0.2. In all, the 25 rejected points contributed 981 to 
the total x 2 , i-e. , an average A\ 2 of about 39 per point! 



For the LHC 
[5] predict 



the In s fits shown in Fig. [4] and Fig. 



a 



VP 

Ppp 



= 107.3 ± 1.2 mb, v / s = 14TeV (22) 
= 0.132 ± 0.001, (23) 



where the quoted errors are due to the uncertainties in 
the fit parameters. 

Deep inelastic scattering (DIS). Berger, Block and Tan 
[3J analyzed the x dependence of the DIS proton structure 
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functions F^x, Q 2 ) by beginning with the assumption 
that the x dependence at extremely small x should man- 
ifest a behavior consistent with saturation of the Froissart 
bound on hadronic total cross sections Q , as is satisfied 
by data on jp, n^p, and pp and pp interactions. They 
treated DIS ep scattering as the hadronic reaction ~f*p, 
demanding that the total cross section o-(j*p) saturate 
the Froissart bound of In 2 s. For the process 7* +p — > X, 
the invariant s is given by s — Q 2 /x, for small x, when 
s 3> m 2 , where m is the proton mass, x is the fractional 
longitudinal momentum of the proton carried by its par- 
ton constituents and Q 2 is the virtuality of the 7*. Thus, 
saturting the Froissart bound [H demands that ff (x, Q 2 ) 
grow no more rapidly than In 2 (1/ x) at very small x. Over 
the ranges of x and Q 2 for which DIS data are avail- 
able, they show that a very good fit to the x dependence 
of ZEUS data [10] is obtained for x < x P = 0.09 and 
Q— ^> m 2 using the expression 



F%{x,Q 2 ) = {l 



v r f p 
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+ B(Q 2 )ln 2 
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X 1 — Xp 
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(24) 



Their fits to DIS data [10| at 24 values of Q 2 cover the 
wide range 0.11 < Q 2 < 1200 GeV 2 . The value x P = 
0.09 is a scaling point such that the curves for all Q 2 
pass through the point x = xp, at which F2(xp,Q 2 ) — 
Fp ~ 0.41, further constraining all of the fits. Figure 
[6] shows that also in ^*p scattering, the Froissart bound 
appears to be saturated. 




FIG. 6: Fits to the proton structure function data, F%(x, Q 2 ) 
vs. x of the form In 2 s, for 13 values of Q 2 . The data are from 
the ZEUS collaboration [13]. The curves show 13 of 28 global 
fits The vertical and horizontal straight lines intersect at 
the scaling point xp — 0.09, (xp) = 0.41. 



The "Aspen" model. The "Aspen" model is a QCD- 



inspired eikonal model [ll| which naturally yields total pp 
and pp cross sections that go as In 2 s at large energies, i.e., 
they automatically saturate the Froissart bound. From 
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FIG. 7: The elastic differential scattering cross section in 
mb/(GeV/c) 2 vs. |t|, in (GeV/c) 2 , using a constrained Aspen 
Model fit (QCD-inspired theory). The solid curve is the predic- 
tion for the reaction pp — > pp at the LHC, at yfs = 14 TeV. The 
dashed curve is the prediction for the reaction pp — > pp at y/s = 1.8 
TeV, at the Tevatron Collider; the data points are from the E710 
experiment. 



the model, one can extract the total cross section ertot, 
the p- value, the elastic cross section cr e i, the nuclear slope 
parameter B (the logarithmic derivative of do~ c \/dt at 
t = 0) and da c \/dt as a function of t, the squared 4- 
momentum transfer. Shown in Fig. [7] is a constrained 
Aspen-model fit [12[ to da e \/dt vs. \t\, for 1.8 TeV , com- 
pared to E710 data at 1.8 TeV; the fit is excellent. Also 
shown is the prediction for the LHC, at 14 TeV. For more 
details af the model and for results for constrained total 
cross sections, etc., see Ref. fl~2| . 

Cosmic ray predictions. There are now available pub- 
lished 13- lg] and preliminary [l7, 18] p-air inelastic pro- 
duction cross sections (c^lf l ail ) that span the enormous pp 
cms (center-of-mass system) energy range 0.1 < y/s < 100 
TeV, reaching energies well above the Large Hadron Col- 
lider (LHC). Further, we expect high statistics results 
from the Pierre Auger Collaboration [l9| in the near fu- 
ture in this ultra-high energy region. Most importantly, 
we now have available very accurate predictions at cos- 



mic ray energies for the total pp cross section, a pp , from 
fits [20j to accelerator data that used adaptive data sift- 
ing algorithms 21 1 and analyticity constraints [22| that 
were not available in the earlier work of Block, Halzen 
and Stanev [23| ■ Here we take advantage of these new pp 
cross section predictions in order to make accurate pre- 
dictions of the cosmic ray p-air total cross sections, to be 
compared with past and future experiments. 

Extracting proton-proton cross sections from pub- 
lished cosmic ray observations of extensive air showers, 
and vice versa, is far from straightforward 24|. By a 
variety of experimental techniques, cosmic ray experi- 
ments map the atmospheric depth at which extensive air 
showers develop and measure the distribution of X ma x, 
the shower maximum, which is sensitive to the inelas- 
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tic p-air cross section cr'^ lr . From the measured X max 
distribution, the experimenters deduce o- p n ^ aiT . We will 



compare published |13Hl5| and recently announced pre- 
liminary values of cr™ 1 ^ with predictions made from u pp , 
using a Glauber model to obtain cr^f^ from a pp . 

^p-air from the A max distribution: Method I. The mea- 
sured shower attenuation length (A m ) is not only sensi- 
tive to the interaction length of the protons in the atmo- 
sphere (Ap_ a j r ) (25|, with 



a — h\ - h 

J1 m njs\p— alr ™ 



]^ = k ^m i (25) 

p— air p— air 

2 



, the proton mass m in g, 

inel 



(with A m and A p _ a j r in gem 

and the inelastic production cross section a p n ^ l air in mb), 
but also depends on the rate at which the energy of the 
primary proton is dissipated into electromagnetic shower 
energy observed in the experiment. The latter effect is 
parameterized in Eq. (|25j) by the parameter k. The value 
of k depends critically on the inclusive particle produc- 
tion cross section and its energy dependence in nucleon 
and meson interactions on the light nuclear target of the 
atmosphere (see Ref. We emphasize that the goal 

of the cosmic ray experiments is cr p _ a i T (or correspond- 
ingly, Ap_ a i r ), whereas in Method I, the measured quan- 
tity is A m . Thus, a significant drawback of Method I is 
that one needs a model of proton-air interactions to com- 
plete the loop between the measured attenuation length 
A m and the cross section oi^fj^, i.e., one needs the value 
of k in Eq. (|25f to compute (T p n ° ail . Shown in Table Q] are 
the widely varying values of k used in the different exper- 
iments. Clearly the large range of fc-values, from 1.15 for 
EASTOP [13 to 1.6 for Fly's Eye differ significantly 
thus making the published values of o- p n ° ail . unreliable. It 
is interesting to note the monotonic decrease over time in 
the fc's used in the different experiments, from 1.6 used 
in Fly's Eye in 1984 to the 1.15 value used in EASTOP in 
2007, showing the time evolution of Monte Carlo models 
of energy dissipation in showers. For comparison, Monte 
Carlo simulations made by Pryke [26| in 2001 of several 
more modern shower models are also shown in Table |U 
Even among modern shower models, the spread is still 
significant: one of our goals is to minimize the impact of 
model dependence on the cr™^ ir determination. 

^p^air fr°m the Amax distribution: Method II. The 
HiRes group [16[ has developed a quasi model-free 
method of measuring o- p n ° aij: directly. They fold into their 
shower development program a randomly generated ex- 
ponential distribution of shower first interaction points, 
and then fit the entire distribution, and not just the trail- 
ing edge, as done in other experiments [TsMlHl ll~7|| . They 
obtain o£^ r = 460 ± 14 (stat) + 39 (syst) - 11 (syst) mb 
at yfs — 77 TeV, which they claim is effectively model- 
independent and hence is an absolute determination [l6| . 

Extraction of a pp from cr p _ ail . The total pp cross sec- 
tion is extracted from ai^f air in two distinct steps. First, 



TABLE I: A table of fc-values, used in experiments and from 
Monte Carlo model simulation 



Experiment 


k 


Fly's Eye 


1.6 


AGASA 


1.5 


Yakutsk 


1.4 


EASTOP 


1.15 


Monte Carlo Results: C.L. Pryke 


Model 


k 


CORSIKA-SIBYLL 


1.15 ±0.05 


MOCCA-SIBYLL 


1.16 ±0.03 


CORSIKA-QGSjet 


1.30 ±0.04 


MOCCA-Internal 


1.32 ±0.03 



one calculates the p-air total cross section, cr p _ a i r , from 
the measured inelastic production cross section using 



incl 

Op— air — Op— air 



0~„ 



a 



q — el 



(26) 



Next, the Glauber method [27j is used to transform the 
measured value of o^ aiT mto tne PP total cross section 
a pp ; all the necessary steps are calculable in the theory. 
In Eq. ([26) the measured cross section for particle pro- 
duction is supplemented with o-^_ ail and cr*~^| r , the elas- 



p— air ' 



tic and quasi-elastic cross section, respectively, as calcu- 
lated by the Glauber theory, to obtain the total cross 
section cr p _ air . The subsequent relation between cr p n J! ail 
and a pp critically involves the nuclear slope parameter 
B, the logarithmic slope of forward elastic pp scattering, 
da pp /dt. A plot of B against cr pp , 5 curves of different 
values of cr"^f^ ir , is shown in Fig. El taking into account 
inelastic screening 28]. The reduction procedure from 
a p n -air to o pp is summarized in Ref. [HI ■ The solid curve 
in Fig. [S] is a plot of B vs. a pp — with B taken from the 
"Aspen" eikonal model and a pp taken from the In 2 s fit 
using analytic amplitudes. The large dot corresponds to 
the value of a pp and B at a/s = 77 TeV, the HiRes energy, 
thus fixing the HiRes predicted value of cr p _ aij . . 

Obtaining a pp from o p n ^ alx . In Fig. El we have plotted 
the values of a pp vs. o- p _ aix that are deduced from the 
intersections of the B-a pp curve with the o~ p _ ail curves in 
Fig. [HI Figure [9] furnishes cosmic ray experimenters with 
an easy method to convert their measured o~ p n ° l ail . to o~ pp , 
and vice versa. 

Determining the k value. In Method I, the extraction 
of Ap_ a i r (or (T^air) from the measurement of A m re- 
quires knowing the parameter k. The measured depth 
A max at which a shower reaches maximum development 
in the atmosphere, which is the basis of the cross sec- 
tion measurement in Ref. [l3} , is a combined measure of 
the depth of the first interaction, which is determined by 
the inelastic cross section, and of the subsequent shower 
development, which has to be corrected for. The model 
dependent rate of shower development and its fluctua- 
tions are the origin of the deviation of k from unity in 
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FIG. 8: B dependence on the pp total cross section a vv . 
The five curves are lines of constant <x™* ir , of 414, 435, 456, 
499 and 542 mb — the central value is the published Fly's Eye 
value, and the others are ±lcr and ±2o\ The solid curve is a 
plot of a QCD-inspired fit of B against a vv , obtained from a 
In 2 s fit. The large dot is the prediction for a l pl^ T at *Js = 77 



TeV, the HiRes energy. 




FIG. 9: A plot of the predicted total pp cross section <7pp, in 
mb vs. the predicted production p-air cross section, Op^ir, 
in mb, where 



Eq. (|25|) . As seen in Table HI its values range from 1.6 for 
a very old model where the inclusive cross section exhib- 
ited Feynman scaling, to 1.15 for modern models with 
large scaling violations. 

Adopting the same strategy that earlier had been used 
by Block, Halzen and Stanev [23$ . we matched the data 
to our prediction of o- l ™ l ail (s), extracting a common value 
for k, neglecting the possibility of a weak energy depen- 
dence of k over the range measured, found to be very 
small in the simulations of Ref . [2] • By combining the 
results of B from the "Aspen" model and Fig.Hl we ob- 
tain our prediction of cr 1 p° l air vs. -v/s, which is shown in 
Fig. 1101 Leaving as a free parameter, we make a x 2 



fit to rescaled o-'^ T (s) values of Flyjs Eye, [HaGASA 
Hi, EASTOP [nl and Yakutsk [lj|, the experiments 
that need a common fc-value. 



Figure [10] is a plot of a 



inel 
p— air 



•■/s, the cms energy 



in GeV, for the two different types of experimental ex- 
traction, using Methods I and II described earlier. Plot- 
ted as published is the HiRes value at yfs = 77 TeV, 
since it is an absolute measurement. We have rescaled in 
Fig.[TO]the_published values of cr"_ e i ir for Fh/s Eye [H, 
AGASA Q, Yakutsk [lj] and EAS-TOP [II], against 
our prediction of cr^fair- using the common value of 
k = 1.264 ± 0.033 ± 0.013 obtained from a v 2 fit, and 
it is the rescaled values that are plotted in Fig. [TQ1 along 
with the rescaled values of ARGO-YB J [lj[ , which were 
not used in the fit. The error in k of 0.033 is the sta- 
tistical error of the x 2 n t- whereas the error of 0.013 is 
the systematic error due to the error in the prediction 



of a. 



inel 
p— air* 



Clearly, we have an excellent fit, with com- 



plete agreement for all experimental points. Our anal- 
ysis gave x 2 — 3-19 for 11 degrees of freedom (the low 
X 2 is likely due to overestimates of experimental errors) . 
We note that our /c-value, k = 1.264 ± 0.033 ± 0.013, 
although somewhat too small for the very low energy 
ARGO-YBJ data, is about halfway between the values 
of CORSIKA-SIBYLL and CORSIKA-QSGSjet found in 
the Pryke simulations [26| , as seen in Table U for model 
predictions. We next compare our measured k parame- 
ter with a recent direct measurement of k by the HiRes 
group 31]- They measured the exponential slope of 
the tail of their X m distribution, A m and compared it 
to the p-air interaction length A p _ a - r that they found. 
Using Eq. (|25|) . they deduced that a preliminary value, 
k = 1.21 + 0.14 — 0.09, in agreement with our value. 




1000 10000 

J5, in GeV 



FIG. 10: A x 2 fit of the renormalized AGASA, EASTOP, 
Fly's Eye and Yakutsk data for cr 1 ™ 1 ^, in mb, as a function of 
the energy, y^s, in GeV. The result of the fit for the parameter 
k in Eq. Jglj) is k = 1.263 ± 0.033. The HiRes point (solid 
diamond), at y^s = 77 GeV, is model- independent and has 
not been renormalized. The renormalized ARGO-YBJ data 
were not used in the fit. 
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Taken together with the goodness-of-fit of our fitted k 
value and the fact that our value, k = 1.264±0.033±0.013 
is compatible with the range of k values from theoretical 
models shown in Table HI the preliminary HiRes value 
31] of k — 1.21 + 0.14 — 0.09 is additional experimen- 
tal confirmation of our overall method for determining 
fc, i.e., our assumptions that the k value is essentially 
energy-independent, as well as being independent of the 
very different experimental techniques for measuring air- 
showers. Our measured k value, k = 1.264±0.033±0.013, 
agrees very well with the preliminary fc-value measured 
by the HiRes group, at the several parts per mil level; in 
turn, they both agree with Monte Carlo model simula- 
tions at the 5-10 part per mil level. 

It should be noted that the preliminary EASTOP mea- 
surement [I?} — at the cms energy = 2 TeV — is at 
an energy essentially identical to the top energy of the 
Tevatron collider, where there is an experimental deter- 



mination of (jpp [32(| , and consequently, no necessity for 



an extrapolation of collider cross sections. Since their 
value of cTp" ^. is in excellent agreement with the pre- 
dicted value of crj^ir) this anchors our fit at its low en- 
ergy end. Correspondingly, at the high end of the cosmic 
ray spectrum, the absolute value of the HiRes experimen- 
tal value of 4 n _!^ r at 77 TeV— which requires no knowl- 
edge of the k parameter — is also in good agreement with 
our prediction, anchoring the fit at the high end. Thus, 
our (T^fair predictions, which span the enormous energy 
range, 0.1 < s/s < 100 TeV, are consistent with cosmic ray 
data, for both magnitude and energy dependence. 

Shown in Fig. [TT] are all of the known pp and pp totl 
cross section data, including the cosmic ray points, span- 
ning the energy region from 2 GeV to 80 TeV, fitted by 
the same In 2 s (saturated Froissart bound) fit of the even 
cross section cr = (°~pp + o'p p )/2. 

CCCR: Cyclotrons to Colliders to Cosmic Rays. We 
have finally reached our goal. This long energy tale of 
accelerator experiments, extending over some 55 years, 
from those using cyclotrons to those using synchrotons 
and then, finally, to those using colliders, has now been 
unified with those experiments using high energy cosmic 
rays as their beams. The accelerator experiments had 
large fluxes and accurate energy measurements, allowing 
for precision measurements: the more precise, the lower 
the energy. On the other hand, the cosmic ray exper- 
iments always suffered from low fluxes of particles and 
poor energy determinations of their events, but made up 
for that by their incredibly high energies. 

The ability to clean up accelerator cross section and 
p-value data by the Sieve algorithm, along with new fit- 
ting techniques using analyticity constraints in the form 
of anchoring high energy cross section fits to the value 
of low energy pp and pp experimental cross sections (and 
their energy derivatives) have furnished us with a preci- 
sion fit — using the a In 2 s form that saturates the Frois- 



a nn even (log (v/m)) 
Fly's Eye 
AGASA 
PP 

pbar-p 
HiRes 




10 2 10 3 
ys, in GeV 



FIG. 11: All known [3flJ a p 
tions, shown together with a vv deduced from the AGASA, Fly's 
Eye and HiRes cosmic ray experiments, pp and pp accelerator to- 
tal cross sections, in mb, vs. the cm. energy ^/s, in GeV. The 
circles are pp and the open squares are pp data. The solid curve is 
a plot of (r o , the even nucleon-nucleon cross section, taken from an 
analytically constrained global In 2 s fit which included the cosmic 
ray data. The AGASA data are the triangles, the Fly's Eye point 
is the diamond and the HiRes point is the open diamond. 



sart Bound — which allows us to make accurate extrapo- 
lations into the LHC and cosmic ray regions, extrapola- 
tions guided by the principles of analyticity and unitarity 
embodied in the Froissart Bound. 
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